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Abstract 

I review what lattice QCD simulations have to tell us about the glueball spectrum. 
We see that the various lattice calculations are in good agreement with each other. 
They predict that prior to mixing with nearby flavour singlet quarkonia the lightest 
glueball states are the scalar at 1.61 it 0.15 GeV, the tensor at 2.26 it 0.22 GeV, and 
the pseudoscalar at 2.19 ± 0.32 GeV. 
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1 Introduction 



My topic here concerns the glueball spectrum. The physics question is: where, in the 
experimentally determined hadron spectrum, are the glueballs hiding? Ideally I should be 
telling you what happens when you simulate QCD with realistically light quarks. But it is 
going to be a few years yet before I can do that. What I will focus on here are lattice glueball 
calculations in the SU(3) gauge theory without quarks. We now know what are the lightest 
states in the continuum (rather than lattice) theory; and I will tell you what they are. The 
next step, if we want to make contact with the real world, is to introduce the physical GeV 
mass scale. This introduces uncertainties which I will try to estimate for you. The final step 
is to discuss possible mixing scenarios with nearby flavour singlet quarkonia. At this stage 
we can look at the experimental spectrum and pinpoint the experimental states most likely 
to have large glueball components. I will not have time to say much about the latter topics 
and refer you instead to ref and ref where you will also find a more complete set of 
references. 

The states in the pure SU(3) gauge theory are glueballs by definition - we only have gluons in 
the theory. If you want hadrons with quarks then you can propagate quarks in this gluonic 
vacuum and tie such propagators together so that the object propagating has the appropriate 
hadronic quantum numbers. That is to say, you calculate hadron masses in the relativistic 
valence quark approximation (the 'quenched approximation') to QCD. The spectrum one 
obtains this way is a remarkably good approximation to the observed light hadron spectrum. 
This is not too surprising: one reason we were able to learn of the existence of quarks in the 
first place is because the low-lying hadrons are in fact well described by a simple valence 
quark picture. 

Of course in this theory, with no vacuum quark loops, we don't have mixing between 
quarkonia and glueballs. There is however reason to believe that this mixing is weak in 
the real world - the Zweig (OZI) rule: hadron decays where the initial quarks all have to 
annihilate are strongly suppressed, e.g.in decays. Such a decay may be thought of as 
quarks —>■ glue — > quarks. Glueball mixing with quarks should therefore be a/OZI sup- 
pressed. As should glueball decays into hadrons composed of quarks. The existence of such 
a suppression is supported by a recent lattice calculation ref 0] . 

The picture we have in mind is therefore as follows. The glueballs will only be mildly affected 
by the presence of light quarks. They will, of course, decay into pions etc. but their decay 
width will be relatively small; and there will be a correspondingly small mass shift. Only 
if there happens to be a flavour singlet quarkonium state close by in mass will things be 
very different, because of the mixing of these nearly degenerate states. In this context we 
expect 'close by' to mean within ~ 100 MeV. So we view the glueballs in the pure SU(3) 
gauge theory as being the 'bare' glueballs of QCD which may mix with nearby quarkonia to 
produce the hadrons that are observed in experiments. All this is an assumption of course, 
albeit well motivated. If true it tells us that the glueballs, whether mixed with quarkonia 
or not, should lie close to the masses they have in the gauge theory. So we now turn to the 
calculation of those masses. 
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2 Glueballs in the SU(3) gauge theory 



In lattice calculations (Euclidean) space-time is discretised to a hypercubic lattice, and the 
volume is made finite and (anti)periodic. So the first step is to be able to calculate masses 
reliably on such a lattice hypertorus. The second is to make sure the volume is large enough. 
Assuming this has been done we obtain the mass spectrum, amj(a), of the discretised theory 
in units of the lattice spacing a. What we actually want is the corresponding spectrum 
of the continuum theory, a = 0. To obtain this we proceed as follows. Theoretically we 
know that in this theory the leading lattice spacing corrections to dimensionless ratios of 
physical masses are 0{o?) where a is the lattice spacing ref [Q. So for small enough a we 
can extrapolate our calculated mass ratios to a = using 

^^^^ = "^^S" = °|+c(am)^ (1) 
am2[a) ^2(0) m2\a = 0) 

where m may be chosen to be mi or m2 or some other physical mass: the difference between 
these choices is clearly higher order in - which we are neglecting. In practice I shall use 
am2 = am = a^/a, where a^cr is the confining string tension as calculated in lattice units, 
and ami = amc will be a glueball mass. 

In Fig |I] I show how the calculated mass ratios, for the lightest scalar glueball, vary with 
a^a. As you see, the behaviour is linear - not surprising given the fact that a^o" is indeed 
small for the values plotted. One fits a straight line and obtains the continuum limit as the 
intercept at a = 0, i.e. at a^cr = 0. We obtain in this way the following lightest three masses 
in the continuum limit: 

""""" 3.65 ±0.11 (2) 



mo++ 
m2++ 
mo-+ 



5.15 ±0.21 (3) 

4.97 ±0.58 (4) 

Although we do not have continuum extrapolations for other masses, the UKQCD lattice 
results strongly suggest that glueballs with other J^*-^ are heavier ref [^. 

The values of the glueball masses that we have used are from refs |^, |^, ^ and for our 
sources of the string tension see refs JI], |[. You may recall that a couple of years ago some 
pubhcity was given to an apparent discrepancy between the GFll and UKQCD predictions 
for the 0"*"^ glueball mass. However, as you can explicitly see in Fig 1, the various calculations 
are entirely in agreement with each other. The discrepancy was an illusion: it arose largely 
from different ways of setting the MeV scale. Such differences should be part of the final 
systematic error on the mass: as below. 

3 Glueballs masses in GeV units 

To introduce MeV units we need the string tension in these units. We can do this by 
calculating the mass of the p or nucleon or ... in the quenched approximation and extrapolate 
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Figure 1: The scalar glueball mass: the GFll values (x) and the rest 
extrapolation to the continuum limit is shown. 
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The best linear 



rrip/y/a oi / ^^a oi ... to the continuum limit as we did for the glueballs. We then set 
nip = 770 MeV or = 930 MeV or ... to obtain a-^cr in MeV units. Because the quenched 
spectrum differs slightly from the real world, these estimates differ slightly. This forms part 
of the error. One can estimate the error ad nauseam, as in refs 0, and this leads to an 
estimate 

= 440 ± 15 ± 35 MeV (5) 
where the first error is statistical and the second is systematic. 

We can now use this value in eqns ^ |^ to express our glueball masses in GeV units. We 
obtain 

= 1.61 ±0.07 ±0.13 GeV (6) 



mo+- 

m2++ = 2.26 ± 0.12 ± 0.18 GeV 
mo-+ = 2.19 ± 0.26 ± 0.18 GeV 



(7) 
(8) 
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These, then, are our best lattice predictions for the hghtest glueballs prior to any mixing 
with nearby quarkonium states. In the case of the O"*""*" the focus is naturally on the /o(1500) 
and any scalar lurking in the /j(1700). The tensor focus is naturally on the /2(1900) and the 
G(2150). With the pseudoscalar things are murky: the lattice calculation has huge errors 
and is far from the obvious i(1490) candidate; but in this case topology is important and 
that is a quantity that is sensitive to light quarks in the vacuum. 
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